Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within similar to 30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O -> OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM/MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with similar to 40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O((aq))(+) species. (C) 2011 American Institute of Physics. doi:10.1063/1.3664746

Glycolaldehyde formation via the dimerization of the formyl radical

Glycolaldehyde, the simplest monosaccharide sugar, has recently been detected in low- and high-mass star-forming cores. Following our previous investigation into glycolaldehyde formation, we now consider a further mechanism for the formation of glycolaldehyde that involves the dimerization of the formyl radical, HCO. Quantum mechanical investigation of the HCO dimerization process upon an ice surface is predicted to be barrierless and therefore fast. In an astrophysical context, we show that this mechanism can be very efficient in star-forming cores. It is limited by the availability of the formyl radical, but models suggest that only very small amounts of CO are required to be converted to HCO to meet the observational constraints

Visualizing the orientational dependence of an intermolecular potential

Scanning probe microscopy can now be used to map the properties of single molecules with intramolecular precision by functionalization of the apex of the scanning probe tip with a single atom or molecule. Here we report on the mapping of the three-dimensional potential between fullerene (C₆₀) molecules in different relative orientations, with sub-Angstrom resolution, using dynamic force microscopy (DFM). We introduce a visualization method which is capable of directly imaging the variation in equilibrium binding energy of different molecular orientations. We model the interaction using both a simple approach based around analytical Lennard–Jones potentials, and with dispersion-force-corrected density functional theory (DFT), and show that the positional variation in the binding energy between the molecules is dominated by the onset of repulsive interactions. Our modelling suggests that variations in the dispersion interaction are masked by repulsive interactions even at displacements significantly larger than the equilibrium intermolecular separation

Colored-noise thermostats \`a la carte

Recently, we have shown how a colored-noise Langevin equation can be used in the context of molecular dynamics as a tool to obtain dynamical trajectories whose properties are tailored to display desired sampling features. In the present paper, after having reviewed some analytical results for the stochastic differential equations forming the basis of our approach, we describe in detail the implementation of the generalized Langevin equation thermostat and the fitting procedure used to obtain optimal parameters. We discuss in detail the simulation of nuclear quantum effects, and demonstrate that, by carefully choosing parameters, one can successfully model strongly anharmonic solids such as neon. For the reader's convenience, a library of thermostat parameters and some demonstrative code can be downloaded from an on-line repository

Carrier-mediated magnetoelectricity in complex oxide heterostructures

While tremendous success has been achieved to date in creating both single phase and composite magnetoelectric materials, the quintessential electric-field control of magnetism remains elusive. In this work, we demonstrate a linear magnetoelectric effect which arises from a novel carrier-mediated mechanism, and is a universal feature of the interface between a dielectric and a spin-polarized metal. Using first-principles density functional calculations, we illustrate this effect at the SrRuO$_3$/SrTiO$_3$ interface and describe its origin. To formally quantify the magnetic response of such an interface to an applied electric field, we introduce and define the concept of spin capacitance. In addition to its magnetoelectric and spin capacitive behavior, the interface displays a spatial coexistence of magnetism and dielectric polarization suggesting a route to a new type of interfacial multiferroic

Modular and predictable assembly of porous organic molecular crystals

Nanoporous molecular frameworks are important in applications such as separation, storage and catalysis. Empirical rules exist for their assembly but it is still challenging to place and segregate functionality in three-dimensional porous solids in a predictable way. Indeed, recent studies of mixed crystalline frameworks suggest a preference for the statistical distribution of functionalities throughout the pores rather than, for example, the functional group localization found in the reactive sites of enzymes. This is a potential limitation for 'one-pot' chemical syntheses of porous frameworks from simple starting materials. An alternative strategy is to prepare porous solids from synthetically preorganized molecular pores. In principle, functional organic pore modules could be covalently prefabricated and then assembled to produce materials with specific properties. However, this vision of mix-and-match assembly is far from being realized, not least because of the challenge in reliably predicting three-dimensional structures for molecular crystals, which lack the strong directional bonding found in networks. Here we show that highly porous crystalline solids can be produced by mixing different organic cage modules that self-assemble by means of chiral recognition. The structures of the resulting materials can be predicted computationally, allowing in silico materials design strategies. The constituent pore modules are synthesized in high yields on gram scales in a one-step reaction. Assembly of the porous co-crystals is as simple as combining the modules in solution and removing the solvent. In some cases, the chiral recognition between modules can be exploited to produce porous organic nanoparticles. We show that the method is valid for four different cage modules and can in principle be generalized in a computationally predictable manner based on a lock-and-key assembly between modules

Ab initio van der Waals interactions in simulations of water alter structure from mainly tetrahedral to high-density-like

The structure of liquid water at ambient conditions is studied in ab initio molecular dynamics simulations using van der Waals (vdW) density-functional theory, i.e. using the new exchange-correlation functionals optPBE-vdW and vdW-DF2. Inclusion of the more isotropic vdW interactions counteracts highly directional hydrogen-bonds, which are enhanced by standard functionals. This brings about a softening of the microscopic structure of water, as seen from the broadening of angular distribution functions and, in particular, from the much lower and broader first peak in the oxygen-oxygen pair-correlation function (PCF), indicating loss of structure in the outer solvation shells. In combination with softer non-local correlation terms, as in the new parameterization of vdW-DF, inclusion of vdW interactions is shown to shift the balance of resulting structures from open tetrahedral to more close-packed. The resulting O-O PCF shows some resemblance with experiment for high-density water (A. K. Soper and M. A. Ricci, Phys. Rev. Lett., 84:2881, 2000), but not directly with experiment for ambient water. However, an O-O PCF consisting of a linear combination of 70% from vdW-DF2 and 30% from experiment on low-density liquid water reproduces near-quantitatively the experimental O-O PCF for ambient water, indicating consistency with a two-liquid model with fluctuations between high- and low-density regions

Using metallic noncontact atomic force microscope tips for imaging insulators and polar molecules: tip characterization and imaging mechanisms

We demonstrate that using metallic tips for noncontact atomic force microscopy (NC-AFM) imaging at relatively large (>0.5 nm) tip-surface separations provides a reliable method for studying molecules on insulating surfaces with chemical resolution and greatly reduces the complexity of interpreting experimental data. The experimental NC-AFM imaging and theoretical simulations were carried out for the NiO(001) surface as well as adsorbed CO and Co-Salen molecules using Cr-coated Si tips. The experimental results and density functional theory calculations confirm that metallic tips possess a permanent electric dipole moment with its positive end oriented toward the sample. By analyzing the experimental data, we could directly determine the dipole moment of the Cr-coated tip. A model representing the metallic tip as a point dipole is described and shown to produce NC-AFM images of individual CO molecules adsorbed onto NiO(001) in good quantitative agreement with experimental results. Finally, we discuss methods for characterizing the structure of metal-coated tips and the application of these tips to imaging dipoles of large adsorbed molecules. © 2014 American Chemical Society